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Abstract 



We introduce a random intersection graph process aimed at modeling sparse evolving 
affiliation networks that admit tunable (power law) degree distribution and assortativity 
and clustering coefficients. We show the asymptotic degree distribution and provide explicit 
asymptotic formulas for assortativity and clustering coefficients. 
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1 Introduction 

Given non-negative weights x = {xi}i>i and y = {yj}j>i, and a nondecreasing positive sequence 
{r(t)}(>i, satisfying lim(_!._|_oo T{t) = +00, let H^^y be the random bipartite graph with biparti- 
tion V = {vi, V2, ■ ■ ■} and W = {wi,W2, • • • }, where edges {wi,Vj} are inserted independently 
and with probabilities 



Here 6 > a > are fixed numbers. Hx^y defines the random intersection graph Gx^y on the 
vertex set V such that any u,v £ V are declared adjacent (denoted u ^ v) whenever they have 
a common neighbor in Hx^y 

Consider, for example, a library where a new item Wi is acquired at time i, and where a new 
user Vj is registered at time j. User Vj picks at random items from a "contemporary literature 
collection" {wi : aT{j) < i < bT{j)} relevant to time j (the interval {i : aT(j) < i < bT{j)} can 
also be considered as the lifetime of the user vj). Every actor Vj and every item Wi is assigned 
weight yj and Xi respectively. These weights model the activity of actors and attractiveness of 
literature items. Now, assume that up to time t the library has acquired items {wi, . . . , w^^(^i)} =: 
^r*(t)) where : N — N is a given nondecreasing function satisfying lim^^+oo T"*(t) = +00. The 
subgraph Hx,y{t) of Hx,y induced by the bipartition Vt = {vi, V2, ■ ■ ■ , vt} and defines the 

random intersection graph Gx,y{t) on the vertex set Vt- vertices u,v £Vt are declared adjacent 
whenever they have a common neighbor in Hx,y{t). The graph Hx,y{t) represents a snapshot 
taken at time t of the "library" records, while the graph Hx^y shows the complete history of the 
"library". Graphs Gx,y{t) and Gx,y represent adjacency relations (between users) observed up 
to time t and during the whole lifetime of the "library", respectively. Assuming, in addition, 
that X and y are realized values of iid sequences X = {Xj}j>i and Y = {Yj}j>i we obtain the 
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random graph Gx,y and the random graph process {Gx,Y{t)}t>i- The parameters of such a 
network model are the probabihty distributions of Xi , Y\ , the functions r, r* and the cutt-offs 
a<b. 

Random intersection graph Gxy is aimed at modehng sparse evolving affiliation networks that 
admit a power law degree distribution and non- vanishing clustering and assortativity coefficients. 
We first observe that choosing inhomogeneous weight sequences x and y one typically obtains 
an inhomogeneous degree sequence of the graph Gx,y ■ vertices with larger weights attract larger 
numbers of neighbours. Consequently, in the case where the probability distributions of Xi and 
Yi have heavy tails, we may expect to obtain a heavy tailed (asymptotic) degree distribution in 
the random graph Gx,Y- Secondly, we observe that if the set W{t) of items selected by a user 
vt is (stochastically) bounded and the lifetimes of two neighbours of vt, say Vs and Vu, intersect, 
then with a non- vanishing probability Vs and Vu share an item from W{t). Consequently, the 
conditional probability a^gu = ^i^s ~ Vu\vs ~ vt,vt ~ u„), called the clustering coefficient, is 
positive and bounded away from zero. In particular, the underlying bipartite graph structure 
serves as a clustering mechanism. 

Let us compare our model with the model of evolving network considered recently by Britton, 
Lindholm and Turova (2011) [6], (see also [5], [18] |19j). In their model vertices are prescribed 
weights, called social indices, and a vertex vt with social index st creates new edges at a rate 
proportional to st- Clearly, both weight sequences {yt}t>i and {st}t>i have the same purpose of 
modeling inhomogeneity of adjacency relations (hence both models possess a power law asymp- 
totic degree distribution). But the model of Britton, Lindholm and Turova (2011) [6| does not 
have the clustering property. We remark, that the role of a bipartite structure in understand- 
ing/explaining clustering properties of some social networks has been discussed in Newman, 
Watts, and Strogatz (2002) [14J. Furthermore, empirically observed clustering properties of 
real affiliation networks have been reproduced with remarkable accuracy by related models of 
random intersection graphs, see [2], [3]. 

In the present paper we only consider the graph Gx,y- We show the asymptotic distribution of 
the degree d{vt) of a vertex vt as time t — )• +oo. We also obtain explicit asymptotic expressions 
for clustering coefficients at\s,m (^s\t,uj '^u\s,ty for s,t, u — )■ +oo such that s < t < u, and for the 
assortativity coefficient (Pearson's correlation coefficient between degrees of adjacent vertices) 

Estd{vs)d{vt) - Bstd{vs)'E,stdy{t) 

fst = ^ • (2j 

^y\'a.rstd{vs)VaLrstd{vt) 

Here E^j denotes the conditional expectation given the event Vs ~ vt and Varsfd(fs) = Estd^(ws) — 
(Estdv{s))'^ . We remark that (empirical) clustering and assortativity coefficients are commonly 
used characteristics of statistical dependence of adjacency relations of real networks. 
Our results are stated in Section 2. Proofs are given in Section 3. 

2 Results 

Degree. We first present our results on the asymptotic degree distribution in Gx,y- We 
obtain a compound probability distribution in the case where r(t) grows linearly in t (clustering 
regime). For T(t) growing faster than linearly in t, we obtain a mixed Poisson asymptotic degree 
distribution. We denote = EXf , and bk = EY/^. 

Theorem 1. Let b > a > 0. Let T{t) = t. Suppose that EX^ < oo and EYi < oo. For t — )• +oo 
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the random variable d{vt) converges in distribution to the random variable 

Ai 

d^ = y^Kj, (3) 

i=i 

where xi,X2,... are independent and identically distributed random variables independent of 
the random variable Ai. They are distributed as follows. For r = 0, 1, 2, . . . , we have 
r + 1 

P(xi = r) = — — P(A2 = r + l) and P(Ai = r) = E e"^' ^, i = l,2. (4) 
EA2 r\ 

Here Ai = 2(6^/2 _ a^/'i)a^Yi and A2 = 2{a-^/^ - b-^/^)biXi. 

The second moment condition EXf < co of Theorem [l] seems to be redundant. 

Theorem 2. Let b > a > and v > I. Let T{t) = t" , t = 1,2 Suppose that EXf < 00 

and EYi < 00. For t — +00 the random variable d{vt) converges in distribution to the random 
variable A3 having the probability distribution 

P(A3 = r) = Ee-^3^, r = 0,1,2,.... (5) 
r! 

Here A3 = 70261^1 and 7 = Av{b^/'^'' - a^/'^''){a-^/'^'' - b-^'"^"). 

Remark 1. The result of Theorem [2] extends to a more general class of increasing nonnegative 
functions r. In particular, assuming that 

t T-^i2t) 

lim — - = 0, sup 7—— < 00, (6) 

t^+<^T{t) t>i r-i(t) 

and that there exists finite limit 

7*=limt-V2 y: 

ar(t)<i<6T(i) a''"0)<i<f'T(j) 

we obtain the convergence in distribution of divt) to A3 defined by ([5]) with A3 = 7*02611^. Here 

denotes the inverse of r (i.e., r(r~^(t)) = t). 
Remark 2. The function r(t) = tint, which grows slower than any power t^ , u > 1, satisfies 
conditions of Remark 1 with 7* = 4(a~^/^ — b~^^'^)(j3^^'^ — a^^^). Furthermore, the functions 
Ti{t) = e^^ * and T2(t) = e*, that grow faster than any power t^, satisfy conditions of Remark 1 
with 7* = 0. 

Clustering. Our next result. Theorem [3| provides explicit asymptotic formulas for clustering 
coefficients. We note that for s < t < n the conditional probabilities Og^tu^ o:t\su and a„|<.t 
are all different and, given < a < 6, mainly depend on the ratios s/t, s/u and t/u. Denote 
PA ■= Pa{s, t, u) = P{vs ~ vt, Vs ~ Vu, vt ~ Vu) the probability that Vg, vt, Vu make up a triangle. 

Theorem 3. Let b > a > 0. Let T{t) = t. Suppose that EXf < 00 and EY^ < 00. Assume that 
s, t, li — 7- +00 so that s < t < u and [an] < [bs\ . We have 

PA = +o(t-2), (7) 

\/stu \\/au x/bsJ 



PA + apfb2t ^{su) ^'^Stisu 

= ^. , ^2u2u!'-U..^~l/2, . +»(^)- (10) 
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Here 



5t\su = ln(n/t) ln(t/s) + ln(u/t) hi{hs/au) + \n{t/ s) \\i{hs/au) + \v?{hs/au), 
5s\ty_ = \n{u/t)\n{hs/au) + \v?{hs/au), 
5u\st = \'a.{t/ s)\a{hs/au) +\T?{hs/au). 

We remark that the condition \au\ < [bs\ of Theorem [3] excludes the trivial case where = 0. 
Indeed, for s < n, the converse inequality \au] > [bs\ means that the lifetimes of Vg and Vu do 
not intersect and, therefore, we have P{vs ~ Vu) = 0. In addition, the inequality [aii] < [bs\ 
implies that positive numbers 6^^^, 5u\st bounded from above by a constant (only 

depending on a and b). 

Assortativity. Let us now consider the sequence of random variables {d{vt)}t>i- We assume 
that r(t) = t. From Theorem [T] we know about the possible limiting distributions for d{vt). 
Moreover, from the fact that Gx,y is sparse we can conclude that, for any given k, the random 
variables d{vt), d{vt+i), . . . , d{vtj^k) are asymptotically independent as i — +oo. An interesting 
question is about the statistical dependence between d{vs) and d{vt) if we know, in addition, 
that vertices Vs and vt are adjacent in Gxy- We assume that s <t and let s,t ^ +oo so that 
bs — at^ +00. Note that the latter condition ensures that the shared lifetime of Vg and vt tends 
to infinity as s,t — ?■ +oo. In this case we obtain that conditional moments 

Bgtdivs) = Bstdivt) + o(l) = + o(l), (11) 
Estd^ivs) = Bstd\vt) + o(l) = 62 + 0(1), 
Estd{vs)d{vt) = (52- A + 0(1), 

are asymptotically constant. Here A = h^^{2h^ + 2h^ + 4(/i6 — /17)) and 

5i = l + h^\h2 + 2/13), 52 = 1 + h^\3h2 + 6/13 + /14 + 6/15 + 4/i6). 

Furthermore, we denote 

hi = 0261, h2 = asbl^, hs = albjb2'y{Vb - \/a), (12) 

/14 = a^bf^"^, /i5 = a2a^b\b2^'^ {Vb - y/a), 

he = alb%^\Vb - h7 = 4blblj\Vb - V^)2, 



and 7 = 2(a — 6 ^/^). A sketch of the derivation of relations ( 11 ) is given in Section 3 below 
Finally, from ([2]) and (11) we obtain that the assortativity coefficient 

A 

is asymptotically constant. 



rst = l- j2 + 0(1) (13) 



We note that each vertex of the graph Gx,y can be identified with the random subset of W, 
consisting of items selected by that vertex, and two vertices are adjacent in Gx,y whenever 
their subsets intersect. Graphs describing such adjacency relations between members of a finite 
family V = {vi, . . . ,Vn} of random subsets of a given finite set W = {wi, . . . , Wm} are called 
random intersection graphs, see |13| . |15j and [9]. Our graph Gx,y is, therefore, a random 
intersection graph evolving in time. One important application of random intersection graphs, 
defined by random subsets of fixed size, is the model of a secure wireless sensor network that 
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uses random predistribution of keys introduced in [8j . Another potential application of random 
intersection graphs is the statistical analysis and modeling of affiliation networks. For example, 
they are useful in explaining clustering properties of the actor network, see H], Finally, we 
mention that asymptotic degree distribution and clustering properties of random intersection 
graphs have been studied in [I], [2], [7|, [TO], [n],[Il], Hi- 
Concluding remarks. We have shown that the random graph Gx,y admits tunable asymptotic 
degree distribution (icluding the power law) and clustering and assortativity coefficients. An 
interesting problem were to study Gx,y and {Gx,Y{t)}t>i in the case where deterministic cutt- 
offs a < 6 in ([T]) are replaced by random cutt-offs Aj < Bj (so that the lifetime [AjT{j), BjT{j)] 
of an actor vj were random). Furthermore, abrupt cutt-offs can be replaced by some smooth 
cutt-off functions. 

3 Proofs 

We first introduce some notation. Then we prove Theorems [T| [2| [3j The proof of Remark 1 goes 
along the lines of the proof of Theorem [2] and is omitted. 

Throughout the proof limits are taken as t — )• +oo, if not stated otherwise. By c we denote 
positive numbers which may only depend on a, b and r. We remark that c may attain different 
values in different places. We say that a sequence of random variables {Ct}t>i converges to zero 
in probability (denoted Ct = op(l)) whenever limsup^ P(|Ct| > e) = for each e > 0. The 
sequence {Ct}t>i is called stochastically bounded (denoted Ct = Op{l)) whenever for each 6 > 
there exists A''^ > such that limsup^ P(|0| > -^<5) < ^■ 
Time intervals 

Tt = {i: aT{t) < i < bT{t)}, T* = {j : ar(j) < i < 6r(j)} (14) 

can be interpreted as lifetimes of the actor vt and attribute Wi respectively. Here and below 
elements of V are called actors, elements of W are called attributes. The oldest and youngest 
actors that may establish a communication link with vt are denoted vt_ and f . Here 

t- = mm{j : Tj n / 0}, t+ = max{j : Tj n / 0}. 

The event "edge {wijVj} is present in Hxy^ is denoted Wi — )• Vj. Introduce random variables 

Ijj = \wi-¥Vj}i = lit) = ^ lij, L = Lt = ^ Uili, 

jeT*\{t} idTt 

hk{I) = Y.Y^j-^l\ ak{I) = Y,Xfi~'''^ ^CN, (15) 
je/ i&i 

Xij = XiYj / ^/ij , Qxyjt) = Xit Xijmm{l,Xij}. (16) 

i&Tt jeT*\{t} 

We remark, that Ui counts all neighbours of Wi in Hxy belonging to the set V \ {fi}, and Lj 
counts all paths of length 2 in Hxy starting from vt- Introduce events 

At = {Xit < 1, i G Tt}, Bt{e) = {Yj < e^j, j G \ {t}}, e > 0. 

By P and E we denote the conditional probability and expectation given X, Y. The conditional 
probability and expectation given Y is denoted Fx and E^. By P^ and Et we denote the 
conditional probability and expectation given Yt. By dry(C)0 'we denote the total variation 



5 



distance between the probability distributions of random variables C and ^. In the case where 
C,C and X,Y are defined on the same probability space, we denote by dTviC^O the total 
variation distance between the conditional distributions of C and ^ given X, Y. 
In the proof we use the following simple fact. For a uniformly bounded sequence of random 
variables {Ct}t>i (i-e., 3 nonrandom h > such that Vt < Q < h almost surely) we have 

Ct = op{l) B(t = o{l). (17) 

In particular, given a sequence of bivariate random vectors {{4>t,'4't)}t>i, defined on the same 
probability space as X,Y, we have 

dTv{(pt,i't) = op{l) dTv{^t,^t) = o{l). (18) 



3.1 Proof of Theorem [T] 

Before the proof we collect auxiliary results. For r(t) := t and Tt, T* defined in (14), we have 

-1/2^(1/2^, +rt- 1/2, 



71 



:= 2(61/2 



E--1/2 -1/2 , I ■-1/2 

72 := 2(a-i/2 - rV^), 



(19) 



where |r|, \r'\ < c. 

Lemma 1. Let t — t- +oo. Assume that EX^ < oo and EYi < oo. We have 

Ve > PiBtie)) = 1 - o(l), 
t-%{[t^,t+]\{t}) = opil), 
P{d(vt) / Lt) = 0(1), 

P(A) = i-o(i), 

QxY{t) = op{l), EQxYit) = o{l). 
For any integers t > a^^{b + b^^) and i G Tt, and any < e < 1 we have 

|Eai(rt) - ai7iti/2| < cait-^'^, \^hi{T* \ {t}) - hi-f2i^''^\ < chir^'^ 



\ {t}) - 6i72i'/'|%(.) < ci"-\eb'^^ + P.Yi\Y,>eH.}) + c6i 
E|ai(rt) - ai7iti/2| < c^''^ 



1/2^^.1/2 



'1/2 



(20) 
(21) 
(22) 
(23) 
(24) 



(25) 
(26) 
(27) 



Proof of Lemma^ Proof of (20). We estimate the probability of the complement event Bt^e) 
using the union bound and Markov's inequality 

P(Bt{e)) < > ^'•^■) = > ^'*-) ^ e-^{t+/t.)EYiI{Y,>eH-} = o{l). 

t-<j<t+ 

Here we estimate t^/t- < c and invoke the bound EYiI|yj;>^| = o(l), for s — t- +c«. 

Proof of (21). Denote 62(t) = i"^Ei<i<t^/- We note that EYi < oo implies b2{t) = op(l). 
The latter bound in combination with the simple inequality t+ /t_ < c implies ( |21[ ) . 
Proof of (23). Let At denote the complement event to At- We have, by the union bound and 
Markov's inequality, 

Pt(At) < PtiXit > 1) < Y,{it)-%^a2 < ca2t-%^. 

i&Tt i&Tt 
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We obtain the bound Pt{Ai) = o(l), which imphes (23), see (17). 

Proof of (22). In view of (17) it suffices to show that Px{d{vt) / Lt) = op(l). We note 
that d{vt) / Lt if and only if S" > 1, where S = Yl' ^ii^i2^iij^i2j ■ Here we denote Yl' = 
E{^„i,}cT, EjGT* nT- Observing that 

Ex\Wjh2j = ^xPiitPi2tPhjPi2j < alY^Yf / {iii2tj) 
we obtain, by Markov's inequahty. 



Px{d{vt) + Lt) = Px{S > 1) < ExS < a^y.Vi ^y/(ni2jr^ 



(28) 



The simple bound E{i^cT, ik ^ ^ ^P^^^^ E'^'(^i^2jr' < cb2i[t-,t+]\{t}). Now, by (Q 
the right-hand side of (28) tends to zero in probability. 

Proof of (24). Denote Xi = max{Xj, 1}, Yj = max{Yj-, 1}, and let Qxrit) denote the sum (16), 
where Ajj is replaced by Xij = XiYj/^/ij. We observe that EX^ < oo and EYi < oo imply 

:= EXlip{Xi) < oo, := EYiip{Yi) < oo, 

for some positive increasing function <^ : [1, +00) — )• [0, +00) satisfying (p{u) — )• +00 as u — )• +00 
(clearly, (p{-) depends on the distributions of Xi and Yi). In addition, we can choose satisfying 
ip{u) < u and ip{su) < (p{s)ip{u), for s,u> 1. From these inequalities one derives the inequality 
min{l, Aij} < ip{Xi)ip{Yj)/ip(y^). The latter inequality implies 



QxYit) < YtQ*xY{t), Q*xYit) ■■= E 



1 



Furthermore, for i ^ Tt and j G T* we have ij > [atji_ =: t^, and t* — +00 as t — +00. Hence 



EQ3,y(t) < E ^ E ^ 



o 



1 



0(1). 



This bound together with the inequalities QxY{t) < QxYit) < ^tQxy(0 shows (24). 
Proof of (|25|) . These inequalities follow from (19). 

Proof of p6[ ). Denote A = |6 — &i72«^''^|, where 6 denotes the sum bi{T^ \ {t}), but with Yj 
replaced by Yj = Yjl^y^Ke^j}^ j {t}. We have 



\hi{T:\{t])-b^^2i'l^\ls, 



< A < Ai + A2 + Ag 



(29) 



where we denote Ai = |6 - E6|, A2 = \Eb - Ebi{T* \ {t})|, A3 = \EbiiT* \ {t}) - ba2i^/^\. 
Next, we evaluate EAi and A2: 



(EAi)2 < E(A?) < 

jeT*\{t} 

A2< r'/'^Y.i 

j&T*\{t} 



j-^EYj^ < £H^\T*\, 



{Yj>s^j} 



< EYlI^y^^,2t_} 



JeT*\{t} 



(30) 
(31) 



In (30) we first apply Cauchy-Schwartz, then use the linearity of variance of an iid sum, and 
finally apply the inequality VarY,- < EY/ < j'^e'^EYj. Invoking ([25|, (|3o|, (|3T| in ([29|) and 
using (IT9|) and |r*| < ci we obtain (p6|). 
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Proof of (|27|. We write E\ai{Tt) - aijit'^/'^\ < EAi + A2, where Ai := \ai{Tt) - Bai{Tt)\ and 
A2 = |Eai(T() — ai7it^/^|, and invoke the inequalities 

(EAi)2 < EA? = ^ j-i(a2 - aj) < cas 

ieTt 



1/2 

and A2 < cai < cOg , see (25). 



□ 



Inequahty (32) below is referred to as LeCam's inequality, see e.g., [T7]. 



Lemma 2. Let S = Ii + I2 + • • • + In be the sum of independent random indicators with 
probabilities P(Ij = 1) = pi. Let A be Poisson random variable with mean pi + ■ ■ ■ + Pn- The 
total variation distance between the distributions Ps of P\ of S and A 



dTviS,A):= sup |P(S G A) -P(A G A)| < 

AC{0,1,2...} 



(32) 



Proof of Theorem^ Before the proof we introduce some notation. Given X,Y, we generate 
independent Poisson random variables 



lli, 6i! £,3i, i 

with conditional mean values 
Er/i = Xit, 



EA 



li 



J^T*\{t} 

{>^ij-pij), EA2i 



i GTt, r = 
E^3i = Xibij2, 

= Xi62^i-^/^, 



1,2,3, 

EA3i = X,5Mi~^/\ 



Here 



52i = bi{T*\{t})-h, 63i = ba2i'^^ - b, b = mm{bi{T*\{t}),bij2i'^^}. 
Finally, we define ^2?. = + ^li) i ^ Tt and introduce random variables 
Lot = '^'niUi, Lrt = ^,miri, r = 1,2,?,. 



(33) 



We assume, in addition, that given X,Y the families of random variables G Tt} and 

{^ri,^ £ Tt,r = 1,2,3,4} are conditionally independent, and that G Tt} is conditionally 

independent of the set of edges of Hx,y that are not incident to vt- 

We are ready to start the proof. In view of (22) the random variables d(vt) and Lt have the 
same asymptotic distribution (if any) . We shall prove that Lt converges in distribution to . In 
the proof we approximate Lt by the random variable L^t, see (34) and (35) below. Afterwards 
we show that L^t converges in distribution to d^,. 

In order to show thatLt and L^t have the same asymptotic distribution (if any) we prove the 
bounds 



dTv{Lt,Lot) = o(l), 
^Lu-L2t\ =o(l). 



dTv{Lot,Lit) = o(l), 
L2t - Lit = op(l). 



(34) 
(35) 
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Here L2t and L^t are marginals of the random vector {L2t,L3t) constructed in (39) below which 
has the property that L2t has the same distribution as L2t and L^t has the same distribution as 

L3t- 



Let us prove the first bound of (34). We shall show below that 

dTviLt,Lot)lA,<t-%^a2iTt). 



(36) 



From the inequality Ea2(7t) = YlieTt''' ^'^2 < ca2 we conclude that a2{Tt) is stochastically 
bounded. Hence t'^'^Yi a2{Tt) = op(l). This bound and (36) combined with (23) imply 



Now the first bound of (34) follows from (18). It remains to prove (36). We denote = 
Yli=lat\ ^i^i + J2i=k+i ViUi and write, by the triangle inequality, 

keTt 

Then we estimate dxv {L'j^-n L'f^) < dTviVk^^k) < {kt)~^Y^^Xl. Here the first inequality follows 
from the properties of the total variation distance. The second inequality follows from Lemma 
[2] and the fact that on the event At we have pkt = ^kt- 



Let us prove the second bound of (34). In view of (17) it suffices to show that dTviLotjL 



op(l). For this purpose we write, by the triangle inequality, 

drviLot, Lit) < ^ dTv{Ll^i,Ll), 

keTt 

where LI := Yli=iat\ Vi'^i + Sl=i+i ^i^u, and estimate 

dTv{Ll^i,Ll) < dTv{ilkUk,rikiik) < P(??fc / 0)dTv{uk,Cik) 
Now, invoking the inequalities 



P(r?fc / 0) = 1 - e-^« < Xk 



and 



dTv{uk,cik) < 



Pkj 



ier-\{t} 



see (32), we obtain from (37), (38) and (24) that 



dTviLo,Li) < Qxvit) = opil). 



(37) 



(38) 



Let us prove the first bound of ( 35 ) . We observe that 

L2t — Lit = 



\L2t — Lit\ 



and 



E m^ii = XI (^^i - i)^{A..>i} ^ QxY{t). 

leTt i&Tt j&T*\{t} 



We obtain B\L2t - Lit\ < EQxvit) = o(l), see (24). 
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Let us prove the second bound of ([35]). We note that the random vector 

(L2t, Lst), L2t = ^ r]i{^4i + A2i), L-st = ^ r/i(^4i + A-si) (39) 

has the marginal distributions of {L2t-,Lzt)- In addition, since A2j, > and at most one of 
them is non-zero, we have |A2i — A3j| = A2i + Aaj. Therefore, we can write 

A := \L2t-ht\ < ^|r/i||A2i- Asil = Vii^2i + Asi) . (40) 

We remark that given X, Y the random variable A2i + Asi has Poisson distribution with (con- 
ditional) mean value 



E(A2i + A3i) = Xii-'/^S,, 6i := \h{T* \ {t}) - 6172^ 



1/2 1 



Therefore, (40) implies EA < f'^^^Yt J2ieTt Next, for < e < 1, we write 



Invoking upper bound (26) for E(5j%^(£) we obtain EIg^(£')A < cb^ a2£ + o{l). Finally, this 
bound combined with Markov's inequality and (20) yields 

P(A > 1) = P({A > 1} n Btie)) + 0(1) < E%(,)A + o(l) < cbl^^a2e + o(l). 

We conclude that P(A / 0) = P(A > 1) = o(l). 

Next we prove that L^t converges in distribution to d^, defined by (|3]). Let 1^ be a random 
copy of Yi, which is independent of X,Y. Given X,Y,Y^, we generate independent Poisson 
random variables rj^, k G Tt with (conditional) mean values E{r]l\X,Y,Y^) = A^*, where Afc* = 
XkY^{kt)~^/'^ . We assume that, given X, y, 1^, the family of random variables {r]^, k € Tt} is 
conditionally independent of {Csk, k G Tt}. Define L* = YlkeT Vk^sk- We note that L* is defined 
in the same way as L^t above, but with Yt replaced by Y^,. Let be defined in the same way 
as d*, but with Ai replaced by A* = y^ai7i. Since L^t has the same distribution as L*, and d^: 
has the same distribution as d^,, it suffices to show that L* converges in distribution to d^. For 
this purpose we show the convergence of Fourier-Stieltjes transforms Ee*^^* — )• Ee*^'^*, for each 
z G (—00, -|-oo). Denote A*[z) = e*^^* — e*^*^*. We shall show below that, for any real z and any 
realized value 1^ there exists a positive constant c* = c*(z, 1^) such that for every < e < 0.5 
we have 

limsup \B{A*{z)\Y^)\ < c*e. (41) 



Clearly, (41 ) implies E(A*(2;)|1^) = o(l). This fact together with the simple inequality |A*(z)| < 
2 yields EA*(z) = o(l), by Lebesgue's dominated convergence theorem. Finally, the identity 
EA*(z) = Ee*^^* - Be"'^* implies Ee*^^* Ee"'^* . 

We fix < e < 0.5 and prove (|41|). Before the proof we introduce some notation. Denote 



Uz) = Ee--\ Uz) = 5]e-^p„ Pr = 'X'^ Y ^^-^te.=r}' ^ = E ^''^^ 

r>0 kdTt kdTt 

6 = iUz) - 1)A - iUz) - 1)A., f{z) = B,e'^^% f{z) = Ee^^^*. 
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Here E denotes the conditional expectation given X, Y, and {^sfc, k G Tt}. By we denote 
the conditional expectation given Y-^. 

Introduce the event V = {\ai{Tt) — 7iait^/^| < et^^"^ min{l, 7iai}} and let V denote the comple- 
ment event. Furthermore, select the number T > 1/e such that P(xi >T)<e. By c^,C2, . . . 
we denote positive numbers which do not depend on t. 

We observe that, given Y^, the conditional distribution of d-^ is the compound Poisson distribution 
with the characteristic function f{z) = e^*^^'^^^^^^\ Similarly, given X,Y,Yi, and {^3^, k G T^}, 
the conditional distribution of L* is the compound Poisson distribution with the characteristic 
function f{z) = e^^f'^^^'^~^\ In the proof of (41) we exploit the convergence A — )• A^, and 
Uz) ^ f^{z). 



Let us prove (41). We write 



E,A*(z) = h + h, h = E,A*{z)Iv, h = E,A*(z)%. 



Here I/2I < 2P^(P) = 2P(P) = o(l). Indeed, the bound P(P) = o(l) follows from m^, by 
Markov's inequality. Next we estimate Ii. Combining the identity E^A*(2:) = E^/(z)(e^^— 1) 
with the inequalities < 1 and — 1| < |'S|el''l, we obtain 

<E^|5|el^lfc<c^E^|,5|Ii,. (42) 

Here we estimated el*^' < e®^* =: c\ using the inequalities 

|5|<2A + 2A„ \ = Y,t-^/^ai{Tt)<2K. 

We remark that the last inequality holds provided that event T) occurs. 

Finally, we show that E^l^llx" ^ + c^^^, + c\\^)£ + o(l). To this aim we write 

5 = (Uz) - 1)(A - A.) + (Uz) - Uz))K 

and estimate \6\ < 2|A — A*| + A^|/>f(z) — fH{z)\. The inequality, which holds on the event V, 
|A — A*| < Yi,e implies E^|A — Aj,|Ix) < c^e with C2 := i^. Next we show that 

E.|/^(z) - Uz)\Iv < (4 + c^)e + 0(1). 

We first split 

Uz) - Uz) = e""{Pr - Pr) =Rl-R2 + R3, 
r>0 

and then estimate separately the terms 

r>T r>T 0<r<T 

Here we denote pr = P(><i = r). The upper bound for R2 follows by the choice of T 

\R2\ <Yp^ = >T)<e. 

r>T 

Next, combining the identities pr = {ai{Tt))^'^ YlkeTt ^~^^'^-^k^{(,3k=r} 
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with the inequality ai{Tt) > t"'^/^ai7i/2, which holds on the event P, we obtain 

I I X — ^ — ^ \ — ^ -^k 2 \ — ^ ^k^3k 

r>T ^ keTt ' keTt 

Now, the identity E^X^^s^ = 02^172 implies E^|i?i|Ix) < c^T^^ < c^e. 

Now we estimate R3. We denote p'^ = ai{Tt){ai^it^/'^)~^pr and observe that the inequahty 
\ai{Tt){ai'yit^^'^)~^ — 1| < e, which holds on the event T>, implies 

I E e^'''(.Pr-p'r)\^v<e Pr<e. 



0<r<T 



0<r<T 



In the last inequality we use the fact that the probabilities {pr}r>o sum up to 1. It follows now 
that 

\Rs\Iv<e+ \Pr-Pr\- (43) 

0<r<T 

Next we estimate 

B^\p'^-Pr\ <E^ \p'^ - E^p; I + I B^p'^ - pr I (44) 
where, by the Cauchy-Schwartz and the linearity of the variance of an iid sum, we have 

(E^Ip; - E^p;|)2 < E^l^; - B^p'f < {ai-iit^'^)~^a2{Tt) < ct~^a2a^^ (45) 



In (45) we first apply the Cauchy-Schwartz inequality, then use the linearity of variance and 
the simple inequality VarXfcI|^g^=r} < 02- In (46) we use the identity E^X/fcI|fg^=^} = aiPr and 
([l9|. From (|44|, ([isf, ^ we conclude that E^|p; - pr\ = ©(t"^/^). Now (|43) implies 



□ 



3.2 Proof of Theorem [2] 

Here we assume that r(t) := . In the proof below we apply the following simple approximations 
^ ^-(i-2.)/(2.) ^ ^1/2^/ ^ Y r'/' = /t^'"^"72 + r'/t-(2-)'\ (47) 

:= 2^.(6(2-)-^ - a(2-)-^), 7^ := 2(a-(2-)-^ - 



where |r'| < c. We also make use of relations (20), (22), (23) and ( |24[ ), which remain valid 
in the case where r(t) = f^, and of the inequalities, for k G Tj, 

|E6i(r,* \ {f}) - 6172^:^/^''^^ I < c6iA;-i/(2.)^ (48) 
m{T*k \ {t}) - bi^'2k'/^'''^\lB,is) < cA:i/(2-)(e6i/2 + EFi V,>,.,_}) + c6ifc-V(2-). (49) 



We note that (|48j) follows from the second identity of (47), and (49) is obtained in the same way 
as (26) above. 
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Proof of Theorem^ Before the proof we introduce some notation. Given e G (0, 1), denote 



Given X, Y, we generate independent Poisson random variables 'i]k,(,3k, k €z Tt, with (conditional) 
mean values E?]^ = A^^, £^3^ = /J^ftiX^ and independent Bernoulli random variables 1^, k £ Tt 
with success probabilities 

p(ifc = 1) = 1 - p(ifc = 0) = Ck. 

We assume that, given X, Y, the sequences {Ifc, k G T^}, {1^, k £ Tj}, {r]k, k £ Tt}, {S,3k, k G Tt} 
are conditionally independent. Next, we introduce random variables 

Lst = ^ ilkisk, Lit = ^ IfcCsfc, L^t = ^ Ifelfc^sfc, Let = ^ Ifclfc- 
fceTt fceTt fceTt fceT* 

Furthermore, we define the random variable L-jt as follows. We first generate X, Y. Then, given 
X,Y, we generate a Poisson random variable with the conditional mean value C,. The realized 
value of the Poisson random variable is denoted Ljt. Thus, we have ^{Ljt = r) = Ee~^C''/r!, 
for r = 0, 1, . . . . 

Now we are ready to prove Theorem [2] In the first step of the proof we show that random 
variables d{vt) and L^t have the same asymptotic distribution (if any). Here we proceed as in 
the proof of ([si]), (|35]) above and make use of ((20]), ([22j), ([23]), ([24]), (|48j), (|49|. In the second 
step we show that L^t converges in distribution to A3. For this purpose we prove that 



dTv{LM,Lu) = 0(1) 
dTv{Lm,Ljt) = 0(1) 



E(L4t - Lst) = 0(1), 



(50) 
(51) 



for every —00 < z < +00, and that there exists c > 0, depending only on a,b^v^ such that for 
any e G (0, 1) we have 

dTv{L5t,LGt) < ca2ble. (52) 



Let us prove (50), (51), (52). The first bound of (50) is obtained in the same way as the first 
bound of ( 34 ) . To show the second bound of ( 50 ) we write 

E(Li5 - Lte) = J] (1 - I'kmktMsk = Ythit-^l^ ^ (1 - tk)Xlhk-^'^ 

k&Tt kGTt 

and apply the simple inequality 

EX|(l-4)<EX,2(l-i;), keTt. (53) 

Here we denote t = mm{k : k £ Tt}. We obtain 

E{Ut - L^t) = EE(L4t - Lst) < StblEXlil - I't) = o(l). 

Here we denote St = T^teTt l^kk^'^^^ and use the simple inequality St < c. Furthermore, 

we invoke the bound EXf(l — I^) = o(l), which holds since t — t- +00 as i — )• +00 

Let us prove (52). Proceeding as in ti7\, tiSv and using the identity Ik = IfcF, we write 



dTviL5t,LQt) < Yl ^'k^i^k + 0)JTy(l3fc,ifc). 
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Next, we estimate I'^drvlCsfc) Ifc) < C|) by LeCam's inequality (32), and invoke the inequality 
P(Ifc / 0) < Xkt- We obtain 

drviL^t, Let) < ^ Ifc^fctCl ^ ^ X/ 

Here we estimated C| < ^Cfc- Now the inequalities 

drviL^t, LQt) < EdTviL^t, Let) < £ ^ EAfctCfc < a2&i5'te 



and St < c imply (52). 

Let us prove the first relation of (51). In view of (17) it suffices to show that dTvi^et, L^t) 
op(l). For this purpose we write 

dTviLet, Ln) < ^AidrviLet, Ljt) + 



where = op(l), see (23), and estimate using LeCam's inequality (j32j) 

lAjrviLeuLn) < IA^ J2 ^'(^'^^fc = 1% < Y^hlt-^ ^ k~^fkXt = op{l) 

Here we used the simple inequality t~^Ylf^^rp^ k^^(3'^X^ < ct~^^ Ylk<bt^ -^t ^"^^ ^^'^^ t^&i 
'EiX'l < oo implies the bound Ylk<n-^k ~ '^^'(1)' as n — )• +oo. 



Finally, we show the second relation of (51). We write Ee*^'^'^' = e^^^'^ and use the bound 

Ytha2^-C = opil). (54) 

We note that, for any real z, the function bounded and uniformly continuous 

for n > 0. Therefore, (54) implies the convergence 



It remains to prove (54). We note that (53) implies 



(55) 



kGTt 



Next, we split 7 = 7x72 and invoke the expression for 7^ obtained from (47). We obtain 

Ytha2^ = Ythiiit-^l^ Yl k-^^~^^'^'^^"^a2 + op{l). (56) 



feGTt 



We observe that (54) follows from ( |55[ ), ( [56^ and the bound 

R := t-V2 J2 (a2 - X^k^^'^r'-' = op{l). 

keTt 



(57) 



In the proof of (57) we use the standard truncation argument. Let e > and let R be defined 
as R above, but with X| replaced by X| = X'^I^x'^<e^k} and 02 replaced by EX|. We have 
R = R + op(l) and P(^ > e^/^) < g-iEi^^ < ce. Letting e ^ we obtain R = op(l). 

□ 
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3.3 Proof of Theorem 3 

Before the proof we state an auxiliary lemma. 

Lemma 3. Denote If = I|x.>ji/2| and = I^y^yj^/^y We have 

A,,(l -If- Ip < min{l, Xi,} < Xij (58) 
Proof of Lemma^ The inequality ^{x-.yiy < If + Ij implies 

Xijil - If - ip < Xij - {Xi, - l)I{A,,>i} = min{l, A,,}. 

□ 

Proof of Theorem]^ The proof of ([9]), ([S]), ((lO]) is very much the same. Therefore, we only prove 
and (H). 

Before the proof we introduce some notation. Denote 



[fit 



T.nTt, Ttu = TtnTu, Tstu = TsnTtnTu, T = TsUTtUTu. 



An attribute Wi is called witness of the edge Vj ~ Vf. whenever lijiik = 1- In this case we say 
that witness Wi realizes the edge vj ~ Vk- Let Ai = {3i : Ij^Iiflj^ = 1} denote the event that all 
three edges of the triangle Vs,vt, Vu are realized by a common witness. Let A2 denote the event 
that all three edges are realized by different witnesses, 

A2 = {3 distinct i,j,k such that lislu = I, Ijslju = 1, Ikthu = ^- 

Let A = {vs ~ vt, Vs ~ Vu, vt ~ Vu} denote the event that vertices Vs,vt, v^ make up a triangle. 
Introduce events Tit = {vg ~ vt,vt ~ Vu} and ICt = {3i / j : liths^jt^ju = 1}, and random 
variables 

S = Iks^kAu: Q = hshthu^js^jt^ju, 

au<k<bs au<i<j<bs 

St = hths^jt^ju, Qt = hsht^jt^ju^kt^ks^rt^ru- 

ik,r)&I,{k,r)f^(i,j) 

Here I denote the set of all ordered pairs & T x T such that i 7^ j. We remark that every 
(i, j) indicates a pair {wi, Wj) of possible witnesses of edges Vg ~ vt and vt ~ Vu respectively. 
We note that for < s < t < u satisfying \au\ < [bs\ the ratios t/s,u/t,u/s G [l,b/a]. Hence 
the variables s, t, u — )• +00 are of the same order of magnitude. 

Let us prove Q. We observe that Ai C A C Ai U A2. Hence 

P(Ai) <P(A) <P(Ai) + P(A2). (59) 

Next, by inclusion exclusion, we write S — Q < Iai ^ S and estimate 

< P(Ai) < E5. (60) 
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Finally, combining ( 59 ) and ( 60 ) with the relations 
BS = 



au<k<bs 



au<i<j<bs 



' stu 



MM 

stu 



au 



'hs 
1 



au<i<j<bs 



^3/2 j 3/2 



- 0{r') 



(61) 



P(A2 



< E 



(62) 



we obtain asymptotic expression ^ for = P(A)- We note that in the first step of (61) we 
apply Lemma [sj and in the last step of (62) we use the inequality ^jg-p^"^ — ^■ 

Let us prove ([s]). We note that ([s]) follows from ^ and the relation 



V{nt) = P(A) + alhlh2-^5t\su + o{t-^' 



' tJsu 



It remains to show (63). From the identity Ht = U /Cj we obtain 

P(7^t) = P(Ai) + P(/Ct) - P(Ai n Kt). 
Next, by inclusion exclusion, we write St — Qt ^ I/Ct ^ St- These inequalities imply 
ESt - EStil - fej - EQtIv, < EI/cJ©, < P{JCt) < ESf 



(63) 



(64) 



(65) 



Here the event = {Yt < et} and e G (0, 1) is non-random. In the remaining part of the proof 
we show that 

1 



BSt = aibib2—=6tis,u + o(i' ), 

t\/SU ' 



(66) 

P(Ain/Ct) = 0(t-3), (67) 
and that there exists c* > which does not depend on s, t, u and e such that, for any e € (0, 1), 
EQilc, < c*et-2 + 0(t-3), BStil - fcj = o(t-2). (68) 



We observe that ( 63 ) follows from ( 64 ) , (|65]) , ( 66 ) and the bounds ([67]) , ( 68 ) . 

Let us prove (66). Since the product pij := PisPuPjtPju is non zero whenever i G Tgt and j G Tf^, 

we have 

E5t = E PiJ = E 5^ Pii- (69) 

It is convenient to split the set : i € T^f, j € Ttu, i / j'} = Ti U • • • U T4 where 

Ti = iT,t \ Tu) X Ttu, T2 = T,t„ X {Ttu \ T^), 

T3 = {(i,j) : e T'stu,* < i}, T4 = {(«,j) : i,j G Tstu, j < i}. 



and write sum (69) in the form 

E5f = Sti + • • • + S't4, 



•Sffc := ^ Pij. 

(*j)6Tfe 



(70) 
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Now (66) follows from (70) and the relations, for 1 < < 4, 

BStk = 'E KsXitXjAju + o{t-^) = alb^j^ Yl i + (71) 



i<fe<4 (ij)eTfe 



In the first step of ( 71 ) we used Lemma [3j 

Let us prove the first bound of (68). We split the collection of vectors(i, j, k,r) 

Q = { (i, J, /c, r) G such that i^j,k^r and (z, j) / (A;, r)| 

into five non intersecting pieces Q = Qi U • • • U Q5, where 

Qi = {{i,j, k,r) ■.i = k}r\Q, Q2 = {{ij, k,r) ■.i = r}nQ, 
Q3 = {{i,j,k,r) : j = k} HQ, = {{i, j, k,r) : j = r} n Q, 

and Q5 = { (i, j, k, r) : all i, j, k, r are distinct } Pi Q, and write 

Denote Q = {(i, j,r) G : all i,j,r are distinct}. Observing that the typical summand of the 
sum Qti is lislit^jt^ju^rt^ru (sincc i = k), we write 

EQtllx)^ < E Y^ Ks^it^jt^juKtKu^Ve 
< 



sV2i3/2 



VieT / 



Here used inequalities < e and X^ieT I — "^^ Similarly, we prove the inequality 

E(5t4lx'e !^ c'et''^. Furthermore, observing that the typical summand of the sum Qt2 is 
lislithvljdjylksht (since i = r), we write 

2 / x 2 

^ 0302 



VieT / VieT / 



stu " ^ 

In the last step we used inequalities Ytt~^Ix>^ < 1 and ^jg^- ^372 < ct~^/^. Hence, EQt2^Vs 
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0{t-^). Similarly, we prove the bound BQ^Iv, = O(r^). Finally, we estimate 



(ij,fc,r)e(Q5 



In the last step we used the inequality Y^^t ^I©^ < e^. 



Let us prove the second bound of ( 68 ) . We have 



2u2 



2, u 



si 



Let us prove (67). The inequalities I^:* ^ 5't) ^Ai < 5" and S < S, where 



5 = ^Ifc and 1^ = Ikshthu, 
imply P(Ai n JCt) = ElAiI/Ct < ^StS. We show that E^*^ = 0(t"^). We split StS = Si + §2, 

ieT jeT\{i} (jJ,fc)eQ 

and estimate 

ESi < E^ ^ XisXitXjtXjuiKu + ^js) = 0{t 

ieTjeT\{i} 

E/S2 < ES'2 < E XisXitXjtXjuXksXku = 0(i (72) 

Here ^2 is defined in the same way as S2, but with replaced by = Iks^ku- 

□ 



3.4 Proof of (11) 



We only give a sketch of the proof. Let s < t satisfy the inequality \at] < [bs\ . An attribute 
Wi is called witness of the edge Vs ~ vt whenever lu^is = 1. The sums 

est = Y ^■^^ ^ Y hshdjs^jt 

i&TsHTt {i,j}CTsnTt 

count witnesses and pairs of witnesses of the edge Vs ^ vt, respectively. We write, by inclusion- 
exclusion, 

est - Qst < < est (73) 

and note that the quadratic term qst is negligibly small. Hence, we approximate 



I 



estil + op(l)), Pivs ~ Vt) = (1 + o(l))Ee,i. 



(74) 
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Given t and i,j £Tt, we denote T*/. = T* \ {t} and introduce random variables 

k&T^t fcer-nT^* t&Tt {i,j}cTt 

We remark that Lt counts pairs {vs ~ vt;Wi), where Wi is a witness of the edge Vg ~ vt in 
Gx,y, for some Vg £ W \ {vt}. In particular, we have d{vt) < Lt. Similarly, Qt counts all 
triples {vs ~ vt]Wi,Wj), where wi and Wj are distinct witnesses of an edge Vg ^ vt- Note that 
a neighbour Vg of vt, which has k witnesses of the edge {vg ~ vt}, contributes 1 to the number 
d{vt) of neighbours of vt- It contributes k to the sum Lt and it contributes (2) to the sum Qt- 
Hence, we always have 

Lt-Qt< d{vt) < Lt. 

We note that the quadratic term Qt is negligibly small and approximate d{vt) = Lt{l + op(l)). 
Combining this approximation with (74) we obtain, for r = 1,2 and u = s,t, 

^std'^ivu) = {Begty^BegtLl + o(l) and Egtd{vg)d{vt) = (Eegtr^BegtLgLt + o(l). (75) 



Next we evaluate expectations in the right-hand sides of ([75j). A straightforward but tedious 
calculation shows that 

Bcgt = e(i + o(i))/ii, 

EegtLg = Ee,tii + o(e) = 6(1 + o(l))(/ii + /i2 + 2/13), 
EcgtLl = EestL2 + o(G) = G(l + o(l))(/ii + 3/12 + 6/13 + /i4 + 6/15 + 4/i6), 
EegtLgLt = 9(1 + o(l))(/ii + 3/12 + 4/i3 + /i4 + 4/15 + 4/17). 



Here we denote B = (st) ' ln{bs/at). We recall that hi are defined in (12) above. Now (11) 
follows from (75). 
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